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Abstract 

The QCD sum rules for moments of production cross section of bb states in e"*" e~ annihilation 
are extremely sensitive to the values of mb and as for moments of large order n. This enables 
one to extract from the existing data on T resonances the values of these parameters with 
a high precision by using a non-relativistic expansion in 1/n. It is found that the sum rules 
fit the data with GeV) = 0.336 ± 0.011 and mb = 4827 ± 7MeV, where the estimate 

of the errors includes the theoretical uncertainty due to subleading in 1/n terms and the 
experimental uncertainty of the e~^ e~ annihilation cross section above the BB threshold. The 
found value of when evolved in two loops up to the Z mass, gives a^^{Mz) = 0.109±0.001. 
The b quark mass mb corresponds to the ‘on shell’ value appropriate for one-loop perturbative 
calculations. 


1 Introduction 


Precision tests of the Standard Model, which are becoming possible due to experimental 
technique, call for a better precision in understanding the parameters of the Model as we see 
it now. In particular the value of the b quark mass enters as a parameter in the predictions for 
the decays of B hadrons, and its understanding is necessary e.g. for a precise determination 
of the weak mixing element Vcb from the data on inclusive decays of B hadrons. The value 
of as determined from a global £t to the data of the LEP experiments at the Z peak: 
af^{Mz) = 0.125 ± 0.005 ± 0.002 (see e.g. in [|l| and also in 0) is argued0 0 to be 
meaningfully higher than what one finds from extrapolation to the Z mass scale of the values 
of the QCD coupling found in a number of analyses at low energies. If confirmed with further 
improvement of experimental accuracy and of the theoretical understanding, this mismatch 
in the values of ag may signal a contribution of a new physics to decays of Z @’ @ . It should 
be noted however, that at present there exists a variety of low-energy estimates of ag, with 
various degree of compatibility with the LEP value, which reflects the present uncertainty 
in the extraction of ag from the low-energy data. (The most recent and extensive review of 
various fits of the value of as is given by Hinchliffe® 0.) 

The purpose of this paper is to present one more determination of ag from the low- 
energy phenomenology, namely from an analysis of the QCD sum rules for the cross section 
of production of the bb hadronic states in e"*" e~ annihilation. Simultaneously the same 
analysis yields a precision determination of the mass parameter irtb. For the states Xg^ the 
QCD sum rules 0 relate the integral moments of the physically measured quantity Rb = 
a{e^ e~ Xg^)/a{e'^ e~ /i"*"/!”) of the form / ds Rb{s )/to theoretically calculable 
derivatives of the vacuum polarization by the vector current (b'j^b). For high enough n the 
moments are saturated by the lowest vector resonances: T’s, and are essentially not sensitive 
to the uncertainty of the cross section in the continuum above the BB threshold. 

It has been noticed long ago0’ 0 that for large n the theoretical calculation of the mo¬ 
ments within the perturbation theory contains as a parameter ag y/n rather than ag, which 
in a dispersive calculation corresponds to dominance of the near-threshold quark-antiquark 
dynamics at typical velocity v ~ 1/y/n. So that ag y/n ^ ag/v is the familiar Coulomb 
parameter. Therefore at large n the Coulomb effects should be explicitly summed up. On 
the other hand, this behavior implies that high moments are very sensitive to the value of ag 
and thus this value can be extracted with high accuracy, even though the experimental input 
is not very precise. Also, by dimension, the n-th moment depends on mf^, which explains 
the high sensitivity at large n to the quark mass. 
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The limiting factor in considering moments with high n is the growth of the relative 
contribution of the non-perturbative terms. The hrst one, appearing in the sum rules for 
heavy quarks, is proportional to the so-called gluon vacuum condensate (a, J i. This 
term grows approximately as relative to the perturbative oneS S . However in the case of 
the sum rules for the bb states its magnitude is still within about 1% at n = 20 (and becomes 
rapidly more important at higher n). In what follows the range of moments from n = 8 to 
n = 20 will be considered, where the non-perturbative contribution can be safely neglected 
and systematically the leading in 1/n approximation will be used, which, in particular, allows 
summation of the Coulomb terms {ag y/n)^- 

An analysis of the sum rules for the bb states along these lines was done quite some time 
ago0 using then available data and a higher range of m from n 25 to n 40, with 
taking into account the hrst non-perturbative term. Those estimates resulted in evaluating 
the mass parameter mb = 4.80 ± O.OSGeV and in an estimate of ag'- dg = 0.30 ± 0.03 at 
a momentum scale of order 1 GeV. Here the analysis is rehned by calculating the effects of 
running of the coupling dg thereby precisely specifying the scale for dg as a function of n. 
Also an uncertainty due to the next term in the 1/n expansion is estimated, which allows to 
consider a lower range of n, where the results are not affected by an uncertainty in the value 
of the gluon condensate. 

One can present several reasons for this method of determining dg being one of the most, 
if not the most, reliable. Firstly, it is fully justihed within the short-distance QCD and 
does not rely on assumptions about local duality, as one has to, when considering rates at 
hxed energy, like, say, in the annihilation of heavy quarkonia into gluons, or even in more 
traditional processes like the total cross section of e~^ e~ annihilation into hadrons at fixed 
energy, or the total hadronic decay rate of the r lepton]^. Secondly, since the parameter 
dgy/n in the considered range of n is of order one, the effects of the coupling are not just 
small corrections but in fact are dominating. Therefore the value of the coupling can be 
determined with a high precision. Thirdly, the moments with sufficiently high n are only 
very weakly sensitive to the poorly known cross section above the BB threshold. Thus 
even quite conservative assumptions about this cross section are sufficient for performing the 
analysis with an acceptable accuracy. 


^The later method of determining as (see e.g. in 0) was recently criticized by Shifmanlll on the basis 


of possible corrections to local duality in the timelike domain. 
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2 Sum rules 


The QCD sum rules^l-J under discussion arise from considering the vacuum polarization op¬ 
erator induced by the electromagnetic current (b'jfj.b) of the b quarks: 


Piif) = 3^/ e”' {0\T(Ux)U0)) |0> 


( 1 ) 


with Qh = —1/3 being the electric charge of b quark. The imaginary part of Pb{s) is related to 
the physically measurable quantity Rb{s): ImPf,(s) = Pf,(s)/127r. Therefore using dispersion 
relation one can write the n-th derivative of Pb{<f) at = 0 in terms of the n-th integral 
moment of Rb{s): 



where the integral runs over all physical values of s where the cross section for production 
of states containing the bb quark pair in e"*" e~ annihilation is non-zero, i.e. it includes the 
T resonances and the continuum above the BB threshold. On the other hand the same 
derivatives as in eq.@ can be calculated by methods of the short-distance QCD. In particular 
in perturbation theory this can be done by the usual calculation of Feynman graphs for Pb{q‘^), 
which in any hnite order in is equivalent to writing those derivatives in terms of the formal 
integrals with the value of Rb{s) given by the perturbation theory in that order in as'- Rl^{s). 
In this way one arrives at the equations 


Rb{s) 


ds 

qH+I 


= R: 


,pt/ 


ds 

sn+l 


-|- non—pert, terms 


( 3 ) 


which express the QCD sum rules in this particular case. The left hand side of the sum rules 
contains a physically measurable quantity, while the right hand side is calculable theoretically 
including the hrst non-perturbative terms0 0 . 


2.1 Non-relativistic limit in the sum rules 


One can notice that at large n the weight function rapidly decreases with energy, so 

that the integral is dominated by the lowest states of bb. For the perturbation theory integral 
in eq. (|^) this implies that at large n the integral is dominated by the states which are close to 
the quark-antiquark threshold sq = 4m^. Defining E as the energy counted from the quark 
threshold, i.e. s = (2m;, -|- P)^, one can write 


Rbis) 


ds 

„n+l 


= Ri 


2 (2mfe -|- E) dE 1 
{2mb + P)2»^+2 ^ {Amiy 


Rb exp 


E 

mb 


-n 


dE 

2mh 


1 + 0 


n 


( 4 ) 


3 











thus effectively replacing the power weight function by an exponential. Moreover, it is clear 
that the integral is determined by the range of energy where \E\/mi, ~ 1/n, i.e. by the 
non-relativistic region near the threshold. If one also recalls the relation of the energy to the 
velocity v of each of the quarks (in the c.m. system): E = one concludes that the 

relevant range of velocity is given by n ~ so that the 1/n expansion of the integrals in 

the sum rules is equivalent to a non-relativistic expansion in 0. 

Furthermore, in the non-relativistic region the spectral density of the electromagnetic 
current can be expressed through the spectral density of the non-relativistic operator 5(r), 
which is convenient to represent in terms of the matrix element (0|(if —of the Green’s 
function (x|(Ff — £^)“^|y) with the non-relativistic Hamiltonian if 0: 


+ F^Im{0|(ff-£)-'|0) . (5) 

V 671 / 

Here is also included the radiative correction, which comes from distances of order l/rub 
and represents a finite renormalization of the electromagnetic current at the threshold. This 
correction will be further discussed and quantified at the end of this section. 

Substituting the expression (|^) into the exponential integral in eq.(^) one finds a relation0 
of the integral moments to the matrix element K{t) = (0| exp(—ff r)|0) of the Euclidean time 
propagation operator exp(—ff r) at r = n/mf,: 


Rb{s) 


ds 

on+l 


= 1 




Stt 


K 


/ n 
\mb 


l + 0{- 
n 


( 6 ) 


ml 

The latter representation in terms of K (r) enables a detailed analysis of the effects of quark 
interaction in the non-relativistic domain both at the perturbative and the non-perturbative 

level. 


2.2 Summing Coulomb effects 


In the lowest order of perturbation theory, i.e. for free quarks the propagation function has 
the familiar form (see e.g. in |T^, and it should also be taken into account that the reduced 
mass in the bb system is mb/2): 


(x|exp(-ffor)|y) = 


mb 
4:71 T 


3/2 


exp 


mb 

4r 


x-y)^ 


so that 


''■•'""(a 


3/2 


(7) 

( 8 ) 
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which reproduces the large n limit of the moments of Rb{s) calculated from the simple loop 
of Fig.l. 

The QCD interaction of quarks via gluons in the leading in ~ 1/n approximation 
reduces to a Coulomb-like potential 

4a., 


V{r) = 


3r 


(9) 


The perturbation theory expansion in this potential generates an expansion for K{t) in 
powers of the parameter as ^/mbf = a^ ^/n. Since we are aiming at the region of n where 
this parameter is of order one, the interaction @ should be taken into account exactly rather 
than perturbatively. In practice this amounts to hnding the propagation function K (r) in the 
Coulomb problem. In terms of Feynman graphs this is equivalent to summing the diagrams of 
Fig.2 in the non-relativistic region. The time-dependent propagator in the Coulomb held can 
be found® by the inverse Laplace transform of the energy-dependent Green’s function 
at negative energy E = —k'^/rrib, which can be derived from a solution of the Schrodinger 
equation in terms of the singular at z = 0 conhuent hypergeometric function U{a, b; z): 


(r| 


^ mb/ 


-1 


m k 


| 0 ) = —e-^^Til-X)U{l-X, 2 , 2kr) = 


m k 


—kr 


—2krt 


1 + t 


dt , 

( 10 ) 

where A = 2mbas/{3k) is the Coulomb parameter for the potential (|^). Setting r = 0 in the 
integral representation in eq. ([I 0 |) gives a divergent integral. However the divergent term does 
not depend on k and thus does not affect the inverse Laplace transform in the variable 
The result for K{t) has the form® K{n/mb) = K^iji/mb) F{'y), where 7 = 2as y/n/Z is the 
Coulomb parameter in the time-dependent representation, and 


F ( 7 ) = 1 - 1 - 2 a /^7 -|- 




exp 


7 


1 + erf 


7 


( 11 ) 


where, as usual, erf(a:) = (2/a/F) Jq exp(—R)dt. One can readily see that the series in 
eq.(p7]) is well converging, so that this expression can be used for practical calculation. The 
function Elj) sums all the terms of the form (a* a/u)^ in the time dependent propagator and 
thus in the sum rules. One can also easily recognize that in the formal limit n ^ 00 the 
function F{'j) is determined by the hrst term in the sum in eq. ([TT|) : F’( 7 ) —>■ exp( 7 ^) 

which is exactly the contribution of the lowest Coulomb bound state|. In the analysis of the 
sum rules for the bb system we will be restricted by the region, where 7 is of order 1. In 

should be noticed however that at finite 7 each p-th term in the sum in eq.(|ll|) receives contribution 
both from the p-th Coulomb bound state and from the continuum above the threshold. 
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connection with this it can be noticed that nnmerically the function F is large in this region: 
F(l) = 49.1243 ..which explains the high sensitivity of the sum rules to the value of ag- 


2.3 Effects of running Og 


In the previous discussion the effects of running of the renormalized coupling were ignored. 
These have to be included in order to be able to relate the moments of Rh to the value of at a 
specihed scale. Naively, one would estimate that at a velocity v the typical momentum flowing 
through the Coulomb gluons in the graphs of Fig.2 is of order rubV, which in the r domain 
translates into the expectation that the normalization for scales as \Jmb/T = mb/\/n. 
However the Coulomb effects substantially modify this estimate in the region where or 
equivalently 7 , is not small. This can be seen from the simple fact that in the formal limit 
7 —>■ cxo the function K{t) is determined by the lowest Coulomb level, which has the intrinsic 
scale given by the Bohr momentum ks = 2asmb/3. Therefore one should expect that the 
proper normalization scale for in the propagation function K{t) is given by ^Jmb/r h^j), 
i.e. K{t) = Ko{t) h^j))) where the function ^(y) describes the Coulomb effects 

on the scale and the overall normalization of ^(y) depends on the renormalization scheme for 
ag. In particular the dominance of the lowest Coulomb level at large 7 implies that in this 
limit ^(y) oc 7 . 

To quantify the effects of running ag we adopt the MS scheme, in which the ag in the 
interaction potential in eq.(y) is0 a^^{K/r) with k = exp(—C — 5/6) and where C = 
0.5772... is the Euler constant. When expressed in terms of normalized at a hxed 

scale /i, the effect of running ag in the potential @ is reduced to a modihcation of the 
potential of the form: 


V(r) = In ^ + 

3r 3r 271 K ’ 


( 12 ) 


where 6 = 9 is the hrst coefficient in the QCD (3 function. (We deal here with the typical 
scales about 1 GeV, so that the appropriate number of active quark flavors in the (3 function is 
Uf = 3.) On the other hand if the function K{t) is written in terms of 7 (p) = 2a^^{fi)y/n/3, 
the same hrst order effect in the running results in the following modihcation for the function 

Fb)- 


F = FiMb) 


dF bag , 

7 ^-:— In 


Ml) mi 


(13) 


dy 27r 

Finally, the modihcation of the Coulomb function ^(y) is clearly the result of the modihcation 
of the potential in eq.(^). To relate the two expressions it is technically convenient to hrst 
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find the modification of the energy-dependent Green’s function by the perturbation of the 
potential as in eq.([I^) and then to hnd the corresponding change of the time-dependent 
propagator by the inverse Laplace transform. For the energy-dependent Green’s function the 
correction has the form: 




-1 




Aas b as m? 
3 27r TT 




(1 -|- f) (1 -|- u) 

tu 


In^ 


K 


(14) 

where each of the integrals runs from zero to inhnity. The integration over r can be easily 
done explicitly, which however is not as simple for the integration over t and u. However the 
integral can be expanded in powers of A and transformed into the r representation term by 
term, which generates an expansion of the modihcation SF{'y) of the function F in powers of 
7 - (In doing this transformation one only needs the transformation rule for a generic term 
l/k^: 1/k^ — r(p/2)).) This procedure works for all terms of the expansion 
of the integral in eq. (|I^) in powers of A, except for the hrst one with A°, which contains a 
manifestly divergent integral. However this term can be easily calculated directly in terms of 
K{t) using the free evolution function in eq.(^. This amounts to calculating the integral 


mt 


dTi 


16 7r^70 ri/2 _ ^^)3/2 


dr exp — 


rribr r 
4ri (r - Ti) 


Inr = 


Sttt 


bVmJ 2 


(15) 

Gollecting all terms and equating the modihcation of the function F{'j) in eq. ([T^) to that 


calculated from the perturbation of the potential in eq.(|^ we hnally hnd the function h{'y) 
in the following form of expansion in powers of 7 : 


liift( 7 ) = ln( 2 K) + + 


(16) 


where F' = dF/d'y and the coefficients Op are given by the integrals 




dp — 


/o ^0 


dy dz 


{^(p/ 2)/2 - ^( 2 ) - In [(1 - y)(l - z)/{l -yz)]} (- hiyzy 

p\T{p/2) {l-yzf 


(17) 


where ' 0 (x) = d In r(a;)/ dx and the integration variables y and are related to those in eq. ([T^) 
by the obvious substitution: y = z = u/{l+u). The integrals (|T7|) for the coefficients 

ttp can be successfully calculated numerically. The convergence of the expansion in eq. ( 0 ) is 
however somewhat slow. In particular, for calculating h{g) with four digit accuracy at 7 = 1 
one has to keep 21 hrst terms in the expansion. (In the numerical analysis, described in the 
next section, were used 30 terms of the expansion.) 
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Thus we find that the normalization point for ag in the parameter 7 in eq. ( 0 ) is itself 
determined by the value of 7 . Therefore in order to take into account the running of the 
coupling constant in one loop one has to determine the value of 7 in eq.([Tl|) by solving the 
equation 

7 = 7 ’ (IS) 


which can be easily done numerically by iterations. 


2.4 Short-distance radiative correction 

Thus discussed radiative effects arise from distances, which are parametrically larger than 
by a factor \/nh{'^), which are the characteristic distances in the non-relativistic Coulomb 
dynamics of quarks at a velocity v ~ 1/y/n. There is however a radiative effect, which comes 
from the distances of order and is associated with correction to the electromagnetic 


vertex. (A discussion of this point in QED can be found in the textbook |^.) To quantify 


this effect let us consider the graphs shown in Fig.3 for the calculation in up to the one-loop 
order of the spatial part of the vector current ( 67 * 6 ) exactly at the bb threshold, i.e. for the 
quark and the antiquark being at rest in the c.m. system. (It is only the spatial part of 
the current which is non vanishing in the c.m. kinematics.) The sum of the graphs can be 
written as 

(b-fib) (^1 + , (19) 

where the one-loop correction term J after integrating over the angular variables can be 
written in the form of an integral over the Euclidean of the virtual gluon: 

f 4 


J = 


r (iSTJ+“’(‘0 * • 


( 20 ) 


where t = /m^ and the integrand contains the infrared-singular term and the regular 

part tn(t), for which the explicit expression is 

4 


wit) = 


(e - 2 f + 2- 3 vm - 


( 21 ) 


The infrared divergence of the singular term corresponds to the singularity of the first 
Coulomb correction oc ag/v at velocity u = 0, and this term is in fact the one accounted 
for in the Coulomb problem calculationHH . Therefore the additional radiative correction is 
described only by the regular part with w{t). The integral with the w{t) from eq.(|21|) can be 
readily done: 

poo 

/ w{t) dt = -8 , ( 22 ) 







which results in the well known expression for the radiative correction to the matrix element 
squared, used in eq.(|^). What is important here is that the integral is determined by 
virtual momenta of order rrib therefore in contrast with the Coulomb terms this radiative 
correction depends on ag normalized at a scale of order A precise specihcation of this 
scale in the MS scheme requires a two loop calculation of the form factor at the threshold, 
which to the best of my knowledge is so far absent. In view of this we have to resort 
to evaluating the scale within the Brodsky - Lepage - Mackenzie (BLM) scheme0, which 
technically amounts to averaging Int with the weight function w{t) determining the one-loop 
correction. For this average we hnd 


1 ^ _ 1 /o°° w{t) \ntdt 3 

2 " 2 w{t) dt ^ 8 ■ 


(23) 


By the BLM prescription the normalization point for is then given by in the V 

scheme (hxed by the Coulomb potential in the momentum space) and by e^^^rrib = 
g-ii/24^^ ^ 0.632 mb in the MS scheme. Therefore we conclude that within the BLM 
method the appropriate value of in the short-distance radiative correction in eq.@ is 
a^(0.632mb). 


3 Numerical analysis 


Summarizing the discussion of the previous section, the large n perturbative formula for 
the moments of Rb{s), including the effects of running coupling and of the 0 {as) radiative 
correction, can be written as 


Rb(s) 


ds 

on+l 


= 1 


16af*(0.632mfe) 

dvr 






rUb 


l + 0{- 
n 


(24) 

where the function ^(y) is given by eq.(^) and its argument is found by solving the equation 
(p!8|) with the function h determined by eqs. ( 0 ) and (O)- Further perturbative corrections 
to the expression (p4D are at least as small as 0(ag) with no additional enhancement in n. 
On the other hand the region of the parameters to be considered is where Og s/n = 0(1). 
Thus the uncertainty due to further radiative corrections is parametrically the same as that 
from the 1/n terms. 

Another class of corrections to eq. (|^) goes beyond perturbation theory and is associated 
with non-perturbative properties of the QCD vacuum. These corrections grow with n, and 
the hrst one, due to the gluon condensate, is well known. In the leading 1/n approximation 
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this correction amountsS to replacing -^( 7 ) in eq.(^) by the expression 


^( 7 ) i-e(7)-4(- 


mi 


72 


( 25 ) 


where the function ^( 7 ) has a rather complicated formS , which however can be approximated 
in the region 7 < 1.5 by simple exponent: .^( 7 ) ~ exp(—O. 87 ). Using this approximation and 
the value 0 of the gluon condensate ~ 0.012 GeU"^, one can estimate that the relative 

magnitude of the non-perturbative term does not exceed about 1% for n < 20. For this 
reason the present analysis is restricted to moments with n < 20 and the non- perturbative 
term is completely ignored. 

On the lower end the range of n is restricted by the applicability of the 1/n expansion. 
Another practical restriction is that the value of Rb is very poorly known experimentally 
above the BB threshold. To keep this uncertainty well suppressed in the moments it is 
desirable to choose larger values of n. 

In the present analysis as the experimental inputs are used the masses and the e"*" e~ widths 
with their experimental errors for the four lowest T resonances, as given by the Particle Data 
Tables®. In connection with this it can be reminded that the value of R{s) for a narrow 
resonance can be written as 


Qtt 

RRe,(s) = 3 r(Bes ^ e+ e“) Mr,, S(s - Mj,,) , 


(26) 


where a contains the renormalization of the electromagnetic a at energy near the mass of 
the resonance: 

- TOmLT 

with Pem{(f) being the full electromagnetic vacuum polarization operator. In the numerical 
analysis here we use the estimate = l.OTa^ at energy in the region of T resonances. This 
correction turns out to be quite essential in the numerical hts. 

As to the cross section above the continuum threshold, i.e. above the T(4S'), a conservative 
assumption is adopted here that the value of Rh{s) is equal to (1 ± 0.5)/3 starting from the 
mass of T(4S'). I believe that the assigned error of 50% well covers the possible uncertainties 
in the contribution of the continuum cross section to the considered moments. 

For a fit by the formula in eq.(24) these input values were used to calculate the “exper¬ 
imental” left hand side of eq.(24) for n = 8, 12, 16 and 20. The choice of these particular 
values is determined, in addition to the discussed restrictions on the range of n, by the require¬ 
ment that the chosen moments are sufficiently representative of the statistically independent 
input data. In other words, if the values of n were chosen in a shorter range with smaller 
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spacing they would essentially represent less statistically independent values than there are 
experimental inputs. Quantitatively this choice is determined by the condition number of the 
covariance matrix of the values of the moments with respect to the dispersion of inputs. 

The actual £t was done in several different ways as regards handling of the uncertainty 
due to the 0{l/n) terms in eg. P^) and the uncertainty in the continuum cross section. The 
effect of the former terms was parametrized as a factor (1 + c/n) in the right hand side of 
eq.(^^ and the value of the coefficient c was either treated as a £t parameter, or hxed at 
c = 0. Then the difference in the £t values of as and mi, in these two calculations gives the 
estimate of the uncertainty due to the unknown 0{l/n) terms. The 50% uncertainty in the 
cross section in the continuum was treated either as an additional statistical error, or the 
extreme cases were taken as hxed. For the all hts were done in terms of GeV). 

The results of the £t, where the uncertainty in the continuum cross section is treated 
as a statistical error are the following. A three parameter £t for as, mi, and c yields: 
af^ilGeV) = 0.336 ± 0.006, mi, = 4827 ± 4 Mel/ and c = -0.59 ± 0.19 with = q.S at 
the minimum corresponding to the central values of the fit parameters. If c is hxed at c = 0 
and a two parameter ht is performed, the minimum of x2 is reached at GeV) = 0.325 

and mi, = 4280 Mel/. However the value of fh® minimum is = 7.2, thus it would 
not be reasonable to ascribe statistical errors to the ht parameters in this case. However, 
the diherence of the central values for as and in the two hts gives the measure of the 
uncertainty due to the 0(l/n) term, which thus leads to the estimate of the errors: 

af^ilGeV) = 0.336 ±0.011 

rufe = 4827±7Mel/. (28) 

A plot, illustrating these two hts, is shown in Fig. 4. 

If the continuum cross section is hxed at the lower extreme bound: Rb{s) = 0.5/3 at 
a/s > M(4S'), the result of a three parameter ht is a^^{lGeV) = 0.338 ± 0.005, mb = 
4828 ± 4 MeV, c = —0.71 ± 0.16 with = 0.1 at the minimum, which values are within the 
range given by eq.(^). However, if c is hxed at c = 0 the two parameter ht gives a statistically 
unacceptable minimal value x^ = 424 at af^^{lGeV) = 0.313 and mb = 4809 MeV. Setting 
the possible continuum cross section at the upper extreme value: Rb{s) = 1.5/3 does not give 
at all a statistically acceptable ht for either three or two parameters. (In the former case 
one hnds x^ = 43 at the minimum at a^^{l GeV) = 0.307 and mb = 4808Mel/, while a two 
parameter ht gives a totally unreasonable minimal value x^ = 800.) As a result we conclude 
that with all the uncertainty in the continuum cross section and in the 0(l/n) terms in the 
sum rules (^^ taken into account the existing data on T resonances are only compatible with 
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values of a* and mb, which he in the range given by eg. (^ 81 ) . 


4 Discussion 


In the derivation of the sum rules in eg. (^) only a one-loop running of a<j is taken into account. 
Therefore it might be possibly argued that the result for ag found in the present analysis is 
not appropriate for a precision two- or three-loop evolution to higher scales, in particular 
for comparison with other estimates at the Z mass scale. To answer to this argument it 
is instructive to evaluate the actual interval of the normalization scale over which the ag 
evolves in eg.(^) when n changes in the considered range, i.e. from n = 8 to n = 20 . 


Using the central fit values for ag and mb in eg.(^), we find that the normalization point 
n{n) = mb /^/n changes from /i(8) = 1.028 GeV down to /i(20) = 0.99 GeV (for reference: 

the corresponding values of the Coulomb parameter 7 are 7 = 0.625 for n = 8 and 7 = 1.007 
at n = 20). Therefore change in In /i is only about 0.04, so that higher loop effects in the (3 


function would be well less than the error for ag in eg.(^). Naturally, a comparison of the 
present result with other estimates of ag at different scales involves evolution of the coupling 
at large intervals of the normalization scale and should be done using the higher loop effects. 
Considering the evolution in two loops with the starting value of a^^{l GeV) in eg.(p^, one 
hnds: 


= 258± 14MeU , 

= 0.254 ±0.006 , 
af^imb) = 0.185 ±0.003 , 

af^iMz) = 0.109 ± 0.001 , (29) 

where is dehned as the position of the singularity in formal numerical solution of the 
two-loop evolution eguation for ag. These derived values of ag are significantly lower than 
the central values of those found from the r decay®: a^^imT-) = 0.33 ± 0.03, and from 
the LEP data0: a^^{Mz) = 0.125 ± 0.005 ± 0.002, although there is a compatibility at 
the level of about 3 ct due to larger errors of the latter determinations. On the other hand 
the value of ag is in a very good agreement with the one determined from the three-gluon 
decay rates of the T resonances, which corresponds to a^^{Mz) = 0.108 ± 0.001 (statistical 
error only) and which is often doubted as being subject to unknown non-perturbative and 
relativistic corrections, as is discussed in p. I believe that the approach used in this paper 
is intrinsically more accurate, than determining the value of ag from small corrections. The 
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reliability of the considered method can be fnrther improved by experimentally measnring 
the cross section of e"*" e~ annihilation above the T(4S') np to abont 12 GeV c.m. energy and 
by a theoretical analysis of the 1/n and terms. 

The valne of the qnark mass in eq.(PE|) reqnires some specihcation, in view of the intrinsic 
nncertainty of the qnark mass in non-pertnrbative QCD0. Apart from snmmation of the 
Conlomb terms, the snm rnles in eq.(^^ contain the QCD radiative effects of the hrst order 
in as- The on shell mass of a heavy qnark is defined in any hnite order of pertnrbation 
theory and enter as a formal parameter in a calcnlation to that order. In this sense the mass 
parameter mb is the on shell mass appropriate for calculations at the one-loop radiative level 
in QCD. Naturally, the derived on shell mass in perturbation theory depends on the order 
in as in which the calculation is done. This is because a perturbative calculation is justified 
inasmuch as it is determined by the short-distance dynamics and thus it is in fact sensitive 
to an off-shell value of the quark mass m(/i) at short distances, which is a combination of 
the on-shell mass, the coupling constant as(/i) and the momentum scale /x corresponding to 
the distance scale involved in the problem. In particular, for a heavy quark and /x -C m the 
hrst-order relation is m(/x) = m — aas{fi) /x , where the constant a is scheme dependent as 
the off-shell mass m{fi) is, unlike the on-shell mass m. Therefore the on-shell mass derived 
from a short-distance calculation should correlate with the value of the coupling as{fi). This 
correlation is conspicuously present in the considered here analysis of the sum rules (^), and 
can be clearly seen on the plot of Fig. 4. Within the described numerical analysis, one can 
hnd that an uncorrelated with ag mass parameter is ml = mb — 0.56 /x for /x = 1 GeV, 

which is determined with a very high statistical accuracy: ml = 4639 ± 2 MeV, which is due 
to the fact that the sum rules in the considered range of n are sensitive to dynamics at 
distances approximately 1 GeV“^. In view of the scheme dependence of the off-shell mass 
the precise implications of this numerical observation are not quite clear .|^ 
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^It is however interesting to note, in this connection, that the numerical value 0.56 of the coefficient a is 
very close to the value a = 16/(97r) « 0.567, which enters the definition of the off-shell mass by the amplitude 
of Thomson scatteringEB. 
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Figure Captions 


Figure 1. The lowest order graph for vacuum polarization by the vector current of b quarks. 


Figure 2. The type of graphs corresponding to summation of all powers of the Goulomb 
parameter ag/v in the vacuum polarisation in the non-relativistic domain near the threshold. 
Dashed lines denote Goulomb gluons. 


Figure 3. The graphs, describing up to the order the matrix element for production of 
quark and antiquark by the vector current exactly at the threshold. 


Figure 4. The contour plot illustrating the £t of the sum rules to the data. (The uncertainty 
in the continuum cross section above the threshold is treated as a statistical error.) The ellipse 
is the one standard deviation contour for and mb (in MeV) for the three parameter fit 
at the optimal value of c: c = —0.59. The heavy dot corresponds to the minimum of in 
the two parameter £t. 
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